Diurnal oscillations of epigenetic modifications are associated with variation in rhythmic expression of homoeologous genes in Brassica napus.

BACKGROUND: Epigenetic modifications that exhibit circadian oscillations also promote circadian oscillations of gene expression. Brassica napus is a heterozygous polyploid species that has undergone distant hybridization and genome doubling events and has a young and distinct species origin. Studies incorporating circadian rhythm analysis of epigenetic modifications can offer new insights into differences in diurnal oscillation behavior among subgenomes and the regulation of diverse expressions of homologous gene rhythms in biological clocks. RESULTS: In this study, we created a high-resolution and multioscillatory gene expression dataset, active histone modification (H3K4me3, H3K9ac), and RNAPII recruitment in Brassica napus. We also conducted the pioneering characterization of the diurnal rhythm of transcription and epigenetic modifications in an allopolyploid species. We compared the evolution of diurnal rhythms between subgenomes and observed that the Cn subgenome had higher diurnal oscillation activity in both transcription and active histone modifications than the An subgenome. Compared to the A subgenome in Brassica rapa, the An subgenome of Brassica napus displayed significant changes in diurnal oscillation characteristics of transcription. Homologous gene pairs exhibited a higher proportion of diurnal oscillation in transcription than subgenome-specific genes, attributed to higher chromatin accessibility and abundance of active epigenetic modification types. We found that the diurnal expression of homologous genes displayed diversity, and the redundancy of the circadian system resulted in extensive changes in the diurnal rhythm characteristics of clock genes after distant hybridization and genome duplication events. Epigenetic modifications influenced the differences in the diurnal rhythm of homologous gene expression, and the diurnal oscillation of homologous gene expression was affected by the combination of multiple histone modifications. CONCLUSIONS: Herein, we presented, for the first time, a characterization of the diurnal rhythm characteristics of gene expression and its epigenetic modifications in an allopolyploid species. Our discoveries shed light on the epigenetic factors responsible for the diurnal oscillation activity imbalance between subgenomes and homologous genes' rhythmic expression differences. The comprehensive time-series dataset we generated for gene expression and epigenetic modifications provides a valuable resource for future investigations into the regulatory mechanisms of protein-coding genes in Brassica napus.

Domains Rearranged Methylase 2 maintains DNA methylation at large DNA hypomethylated shores and long-range chromatin interactions in rice.

DNA methylation plays an important role in gene regulation and genomic stability. However, large DNA hypomethylated regions known as DNA methylation valleys (DMVs) or canyons have also been suggested to serve unique regulatory functions, largely unknown in rice (Oryza sativa). Here, we describe the DMVs in rice seedlings, which were highly enriched with developmental and transcription regulatory genes. Further detailed analysis indicated that grand DMVs (gDMVs) might be derived from nuclear integrants of organelle DNA (NORGs). Furthermore, Domains Rearranged Methylase 2 (OsDRM2) maintained DNA methylation at short DMV (sDMV) shores. Epigenetic maps indicated that sDMVs were marked with H3K4me3 and/or H3K27me3, although the loss of DNA methylation had a negligible effect on histone modification within these regions. In addition, we constructed H3K27me3-associated interaction maps for homozygous T-DNA insertion mutant of the gene (osdrm2) and wild type (WT). From a global perspective, most (90%) compartments were stable between osdrm2 and WT plants. At a high resolution, we observed a dramatic loss of long-range chromatin loops in osdrm2, which suffered an extensive loss of non-CG (CHG and CHH, H = A, T, or C) methylation. From another viewpoint, the loss of non-CG methylation at sDMV shores in osdrm2 could disrupt H3K27me3-mediated chromatin interaction networks. Overall, our results demonstrated that DMVs are a key genomic feature in rice and are precisely regulated by epigenetic modifications, including DNA methylation and histone modifications. OsDRM2 maintained DNA methylation at sDMV shores, while OsDRM2 deficiency strongly affected three-dimensional (3D) genome architectures.

Integrated 3D genome, epigenome and transcriptome analyses reveal transcriptional coordination of circadian rhythm in rice.

Photoperiods integrate with the circadian clock to coordinate gene expression rhythms and thus ensure plant fitness to the environment. Genome-wide characterization and comparison of rhythmic genes under different light conditions revealed delayed phase under constant darkness (DD) and reduced amplitude under constant light (LL) in rice. Interestingly, ChIP-seq and RNA-seq profiling of rhythmic genes exhibit synchronous circadian oscillation in H3K9ac modifications at their loci and long non-coding RNAs (lncRNAs) expression at proximal loci. To investigate how gene expression rhythm is regulated in rice, we profiled the open chromatin regions and transcription factor (TF) footprints by time-series ATAC-seq. Although open chromatin regions did not show circadian change, a significant number of TFs were identified to rhythmically associate with chromatin and drive gene expression in a time-dependent manner. Further transcriptional regulatory networks mapping uncovered significant correlation between core clock genes and transcription factors involved in light/temperature signaling. In situ Hi-C of ZT8-specific expressed genes displayed highly connected chromatin association at the same time, whereas this ZT8 chromatin connection network dissociates at ZT20, suggesting the circadian control of gene expression by dynamic spatial chromatin conformation. These findings together implicate the existence of a synchronization mechanism between circadian H3K9ac modifications, chromatin association of TF and gene expression, and provides insights into circadian dynamics of spatial chromatin conformation that associate with gene expression rhythms.

The dynamics of lncRNAs transcription in interspecific F1 allotriploid hybrids between Brassica species.

Interspecific hybridization is the intrinsic forces behind genome evolution. Long non-coding RNAs (lncRNAs) are important for plant biological processes regulation. However, it is unclear that these non-coding fractions are impacted by interspecific hybridization. Here we examined the profiles of lncRNAs by comparing them with coding genes in Brassica napus, three accessions of Brassica rapa, and their F1 hybrids. 6206 high-confidential lncRNAs were identified in F 1 hybrids and their parentals, and the lncRNAs transcriptome in the F1 hybrids was reprogrammed by the genome shock. Notably, genome-wide unbalanced of lncRNAs were observed between An and Ar subgenomes, ELD (Expression Level Dominance) was biased toward the An -genome in F1 hybrids, and ELD of non-conserved lncRNAs was more than conserved lncRNAs. Our findings demonstrate that the reprogramed lncRNAs acts as important role in enhancing plant plasticity, leading to the acquisition of desirable traits in polyploid Brassica species.

Circ_0058792 regulates osteogenic differentiation through miR-181a-5p/Smad7 axis in steroid-induced osteonecrosis of the femoral head.

Osteonecrosis of the femoral head (ONFH) caused by steroids is a severe orthopedic disorder resulting from the use of high-dose steroid drugs, characterized by structural changes in the bone, joint dysfunction, and femoral head collapse. CircRNAs and miRNAs have increasingly been suggested to play pivotal roles in osteogenic differentiation and osteogenesis. Significant upregulation of circ_0058792 was observed in patients with steroid-induced ONFH. Bioinformatic analysis showed that circ_0058792 might act as a sponge for miR-181a-5p. This study further investigated the mechanisms underlying the role of circ_0058792 and miR-181a-5p in osteogenic differentiation in methylprednisolone-induced ONFH rats and MC3T3-E1 cells. The results showed a notable decrease in the serum of miR-181a-5p in methylprednisolone-induced ONFH rats. Silencing of circ_0058792 using siRNAs and overexpression of miR-181a-5p significantly increased alkaline phosphatase activity and matrix mineralization capacity. Additionally, markers for osteogenic differentiation were significantly upregulated in miR-181a-5p-transfected cells. However, overexpression of circ_0058792 and the addition of the miR-181a-5p inhibitor reversed this increase. Smad7 was identified to be miR-181a-5p's direct target and circ_0058792 was confirmed to be miR-181a-5p's competing endogenous RNA (ceRNA). Upregulation of miR-181a-5p promotes phosphorylation of Smad2 and Smad3. Furthermore, circ_0058792 and miR-181a-5p had opposing effects on Smad7 expression. Collectively, these findings indicate that circ_0058792 regulates osteogenic differentiation by sponging miR-181a-5p via the TGF-ß/Smad7 pathway. These findings elucidated the functions of circ_0058792 and miR-181a-5p in the regulation of steroid-induced ONFH. Our findings also indicated that circ_0058792 and miR-181a-5p are possible diagnostic markers and therapeutic targets for treating steroid-induced ONFH.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species.

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

Effect of particle size on the oral absorption of isoliquiritigenin nanocrystals

Abstract As one of the most promising formulations for poorly water-soluble drugs, nanocrystals have been attracting increasing attention in recent years. Isoliquiritigenin (ISL) is a flavonoid with a chalcone structure, and possesses many biological activities. However, its clinical application is significantly limited mainly due to its low oral bioavailability caused by poor hydrophilicity. To address this, ISL nanocrystals were developed in this study to improve its oral bioavailability. Three types of nanocrystals with differing particle size; R1, R2, and R3, were prepared by anti- solvent precipitation or anti-solvent precipitation combined with sonication, which was optimized by single-factor experiments. These nanocrystals were characterized based on their physical properties, in vitro release, and in vivo absorption performance. The mean particle size of R1, R2, and R3 was 555.7, 271.0, and 46.2, respectively. The dissolution ratio of ISL in the nanocrystals was significantly improved, with the quickest rate recorded in R2. Peak concentration and area under the concentration-time curve of R2 after oral administration in rats was 5.83- and 2.72-fold higher than that of the ISL solution, respectively. These findings indicate that the dissolution and absorption of ISL can be significantly enhanced by nanocrystals, and the dissolution behavior and pharmacokinetic properties of nanocrystals is significantly influenced by particle size.

Traveling Waves and Estimation of Minimal Wave Speed for a Diffusive Influenza Model with Multiple Strains.

Antiviral treatment remains one of the key pharmacological interventions against influenza pandemic. However, widespread use of antiviral drugs brings with it the danger of drug resistance evolution. To assess the risk of the emergence and diffusion of resistance, in this paper, we develop a diffusive influenza model where influenza infection involves both drug-sensitive and drug-resistant strains. We first analyze its corresponding reaction model, whose reproduction numbers and equilibria are derived. The results show that the sensitive strains can be eliminated by treatment. Then, we establish the existence of the three kinds of traveling waves starting from the disease-free equilibrium, i.e., semi-traveling waves, strong traveling waves and persistent traveling waves, from which we can get some useful information (such as whether influenza will spread, asymptotic speed of propagation, the final state of the wavefront). On the other hand, we discuss three situations in which semi-traveling waves do not exist. When the control reproduction number [Formula: see text] is larger than 1, the conditions for the existence and nonexistence of traveling waves are determined completely by the reproduction numbers [Formula: see text], [Formula: see text] and the wave speed c. Meanwhile, we give an interval estimation of minimal wave speed for influenza transmission, which has important guiding significance for the control of influenza in reality. Our findings demonstrate that the control of influenza depends not only on the rates of resistance emergence and transmission during treatment, but also on the diffusion rates of influenza strains, which have been overlooked in previous modeling studies. This suggests that antiviral treatment should be implemented appropriately, and infected individuals (especially with the resistant strain) should be tested and controlled effectively. Finally, we outline some future directions that deserve further investigation.

Transcriptome Analysis Reveals the Effect of Long Intergenic Noncoding RNAs on Pig Muscle Growth and Fat Deposition.

Muscle growth and fat deposition are the two important biological processes in the development of pigs which are closely related to the pig production performance. Long intergenic noncoding RNAs (lincRNAs), with lack of coding potential and the length of at least 200nt, have been extensively studied to play important roles in many biological processes. However, the importance and molecular regulation mechanism of lincRNAs in the process of muscle growth and fat deposition in pigs are still to be further studied comprehensively. In our study, we used the data, including liver, abdominal fat, and longissimus dorsi muscle of 240 days' age of two F2 full-sib female individuals from the white Duroc and Erhualian crossbreed, to identify 581 putative lincRNAs associated with pig muscle growth and fat deposition. The 581 putative lincRNAs shared many common features with other mammalian lincRNAs, such as fewer exons, lower expression levels, and shorter transcript lengths. Cross-tissue comparisons showed that many transcripts were tissue-specific and were involved in the important biological processes in their corresponding tissues. Gene ontology and pathway analysis revealed that many potential target genes (PTGs) of putative lincRNAs were involved in pig muscle growth and fat deposition-related processes, including muscle cell proliferation, lipid metabolism, and fatty acid degradation. In Quantitative Trait Locus (QTLs) analysis, some PTGs were screened from putative lincRNAs, MRPL12 is associated with muscle growth, GCGR and SLC25A10 were associated with fat deposition, and PPP3CA, DPYD, and FGGY were related not only to muscle growth but also to fat deposition. Therefore, it implied that these lincRNAs might participate in the biological processes related to muscle growth or fat deposition through homeostatic regulation of PTGs, but the detailed molecular regulatory mechanisms still needed to be further explored. This study lays the molecular foundation for the in-depth study of the role of lincRNAs in the pig muscle growth and fat deposition and further provides the new molecular markers for understanding the complex biological mechanisms of pig muscle growth and fat deposition.

Transcriptome Analysis Suggests the Roles of Long Intergenic Non-coding RNAs in the Growth Performance of Weaned Piglets.

Long intergenic non-coding RNAs (lincRNAs) have been considered to play a key regulatory role in various biological processes. An increasing number of studies have utilized transcriptome analysis to obtain lincRNAs with functions related to cancer, but lincRNAs affecting growth rates in weaned piglets are rarely described. Although lincRNAs have been systematically identified in various mouse tissues and cell lines, studies of lincRNA in pigs remain rare. Therefore, identifying and characterizing novel lincRNAs affecting the growth performance of weaned piglets is of great importance. Here, we reconstructed 101,988 lincRNA transcripts and identified 1,078 lincRNAs in two groups of longissimus dorsi muscle (LDM) and subcutaneous fat (SF) based on published RNA-seq datasets. These lincRNAs exhibit typical characteristics, such as shorter lengths and lower expression relative to protein-encoding genes. Gene ontology analysis revealed that some lincRNAs could be involved in weaned piglet related processes, such as insulin resistance and the AMPK signaling pathway. We also compared the positional relationship between differentially expressed lincRNAs (DELs) and quantitative trait loci (QTL) and found that some of DELs may play an important role in piglet growth and development. Our work details part of the lincRNAs that may affect the growth performance of weaned piglets and promotes future studies of lincRNAs for molecular-assisted development in weaned piglets.

Identification and Functional Prediction of Long Intergenic Non-coding RNAs Related to Subcutaneous Adipose Development in Pigs.

An increasing number of studies have shown that long intergenic non-coding RNAs (lincRNAs) are a very important class of non-coding RNAs that plays a vital role in many biological processes. Adipose tissue is an important place for storing energy, but few studies on lincRNAs were related to pig subcutaneous fat development. Here, we used published RNA-seq data from subcutaneous adipose tissue of Italian Large White pigs and identified 252 putative lincRNAs, wherein 34 were unannotated. These lincRNAs had relatively shorter length, lower number of exons, and lower expression level compared with protein-coding transcripts. Gene ontology and pathway analysis indicated that the adjacent genes of lincRNAs were involved in lipid metabolism. In addition, differentially expressed lincRNAs (DELs) between low and high backfat thickness pigs were identified. Through the detection of quantitative trait locus (QTL), DELs were mainly located in QTLs related to adipose development. Based on the expression correlation of DEL genes and their differentially expressed potential target genes, we constructed a co-expression network and a potential pathway of DEL's effect on lipid metabolism. Our study identified and analyzed lincRNAs in subcutaneous adipose tissue, and results suggested that lincRNAs may be involved in the regulation of subcutaneous fat development. Our findings provided new insights into the biological function of porcine lincRNAs.

Time course analysis based on gene expression profile and identification of target molecules for colorectal cancer.

BACKGROUND: The study aimed to investigate the expression changes of genes in colorectal cancer (CRC) and screen the potential molecular targets. METHODS: The GSE37178 of mRNA expression profile including the CRC samples extracted by surgical resection and the paired normal samples was downloaded from Gene Expression Omnibus database. The genes whose expressions were changed at four different time points were screened and clustered using Mfuzz package. Then DAVID was used to perform the functional and pathway enrichment analysis for genes in different clusters. The protein-protein interaction (PPI) networks were constructed for genes in the clusters according to the STRING database. Furthermore, the related-transcription factors (TFs) and microRNAs (miRNAs) were obtained based on the resources in databases and then were combined with the PPI networks in each cluster to construct the integrated network containing genes, TFs and miRNAs. RESULTS: As a result, 314 genes were clustered into four groups. Genes in cluster 1 and cluster 2 showed a decreasing trend, while genes in cluster 3 and cluster 4 presented an increasing trend. Then 18 TFs (e.g., TCF4, MEF2C and FOS) and 18 miRNAs (e.g., miR-382, miR-217, miR-1184, miR-326 and miR-330-5p) were identified and three integrated networks for cluster 1, 3, and 4 were constructed. CONCLUSIONS: The results implied that expression of PITX2, VSNL1, TCF4, MEF2C and FOS are time-related and associated with CRC development, accompanied by several miRNAs including miR-382, miR-217, miR-21, miR-1184, miR-326 and miR-330-5p. All of them might be used as potential diagnostic or therapeutic target molecules for CRC.

Prediction of feature genes in trauma patients with the TNF rs1800629 A allele using support vector machine.

BACKGROUND: Tumor necrosis factor (TNF)-α variant is closely linked to sepsis syndrome and mortality after severe trauma. We aimed to identify feature genes associated with the TNF rs1800629 A allele in trauma patients and help to direct them toward alternative successful treatment. METHODS: In this study, we used 58 sets of gene expression data from Gene Expression Omnibus to predict the feature genes associated with the TNF rs1800629 A allele in trauma patients. We applied support vector machine (SVM) classifier model for classification prediction combining with leave-one-out cross validation method. Functional annotation of feature genes was carried out to study the biological function using database for annotation, visualization, and integrated discovery (DAVID). RESULTS: A total of 133 feature genes were screened out and was well differentiated in the training set (14 patients with variant, 15 with wild type). Moreover, SVM classifier peaked in predictive accuracy with 100% correct rate in training set and 86.2% in testing set. Interestingly, functional annotation showed that feature genes, such as HMOX1 (heme oxygenase (decycling) 1) and RPS7 (ribosomal protein S7) were mainly enriched in terms of cell proliferation and ribosome. CONCLUSION: HMOX1 and RPS7 may be key feature genes associated with the TNF rs1800629 A allele and may play a crucial role in the inflammatory response in trauma patients. Moreover, the cell proliferation and ribosome pathway may contribute to the progression of severe trauma.

Decreased expression of microRNA-130a correlates with TNF-α in the development of osteoarthritis.

OBJECTIVE: Increased expression of tumor necrosis factor a (TNF-α) has emerged as an important inflammatory factor in osteoarthritis (OA) and other joint diseases. The study was performed to investigate whether the expression of TNF-α in human chondrocytes was regulated by miRNAs. METHODS: MiRNA-130a and TNF-α expression in cartilage specimens was examined in patients with knee osteoarthritis, chondrocytes and osteoarthritis rat model. Chondrocytes were transfected with siRNAs as a gene silencing methods. Expression of genes and proteins were analyzed by real-time PCR and western blotting respectively. RESULTS: Increased TNF-α and decreased miRNA-130a were observed in tissues from osteoarthritis patients. Moreover, we found a highly negitive correlation between miRNA-130a and TNF-α. Next, miRNA-130a loss-of-function increased the expression of TNF-α and promoted inflammation in chondrocytes. It was reasonable that miRNA-130a regulated a distinct underlying molecular and pathogenic mechanism of OA by forming a negative feedback loop with TNF-α. Furthermore, there were the abnormalities of bone metabolism in OA rat, which showed the miRNA-130a and TNF-α dysfunction that was one of important factors for the occurrence and development of OA. CONCLUSIONS: Our results indicated that miR-130a played an important role in regulating the expression of TNF-α in human chondrocytes and identified miR-130a as a novel therapeutic target in OA.

Identification of key genes associated with colorectal cancer based on the transcriptional network.

Colorectal cancer (CRC) is among the most lethal human cancers, but the mechanism of the cancer is still unclear enough. We aimed to explore the key genes in CRC progression. The gene expression profile (GSE4183) of CRC was obtained from Gene Expression Omnibus database which included 8 normal samples, 15 adenoma samples, 15 CRC samples and 15 inflammatory bowel disease (IBD) samples. Thereinto, 8 normal, 15 adenoma, and 15 CRC samples were chosen for our research. The differentially expressed genes (DEGs) in normal vs. adenoma, normal vs. CRC, and adenoma vs. CRC, were identified using the Wilcoxon test method in R respectively. The interactive network of DEGs was constructed to select the significant modules using the Pearson's correlation. Meanwhile, transcriptional network of DEGs was also constructed using the g: Profiler. Totally, 2,741 DEGs in normal vs. adenoma, 1,484 DEGs in normal vs. CRC, and 396 DEGs in adenoma vs. CRC were identified. Moreover, function analysis of DEGs in each group showed FcR-mediated phagocytosis pathway in module 1, cardiac muscle contraction pathway in module 6, and Jak-STAT signaling pathway in module 19 were also enriched. Furthermore, MZF1 and AP2 were the transcription factor in module 6, with the target SP1, while SP1 was also a transcription in module 20. DEGs like NCF1, AKT, SP1, AP2, MZF1, and TPM might be used as specific biomarkers in CRC development. Therapy targeting on the functions of these key genes might provide novel perspective for CRC treatment.

Comparison of Clinical Outcomes of Phalangeal Fracture Treated with Dorsolateral Approach or Post-middle Approach Using AO Mini Titanium Plate.

The aim of this study is to investigate the clinical outcomes of various fixation methods for proximal phalangeal fractures with Arbeitsgemeinschaft für Osteosynthesefragen (AO) mini titanium plate by dorsolateral approach or post-middle approach. Clinical results of 62 fingers of 53 patients with proximal phalangeal fracture were evaluated. For dorsolateral approach, the lateral bundle of the extensor tendon was drawn away to expose the fracture part of the bone. After reduction, the plate was located at the dorsolateral side of the bone. For post-middle approach, the extensor tendon was split to expose the fracture part of the bone. After reduction, the plate was fixed to the proximal phalangeal side of the bone, and the extensor tendon was repaired with 3-0 nonabsorbable silk sutures. We found low overall complication rates in both groups. The mean total active motion (TAM) for the dorsolateral group and post-middle group was 234.60° ± 22.63° and 221.08° ± 25.69°, respectively. There was a statistical significance between the two groups (P = 0.037 < 0.05), indicating that TAM was notably affected by various fixation methods. With AO mini titanium plate, movement in dorsolateral approach group was significantly higher than in post-middle approach group. Dorsolateral approach is an acceptable technique of incision for proximal phalangeal fractures.

[Summary of Hui prescriptions for stroke].

Current Hui prescriptions are mostly recorded in the Arabic language. Their fussy and inconsistent names (Arabic names) result in the restriction in the clinical application of Hui prescriptions. Having collected and screened out 101 Hui prescriptions for stroke, the author further studied some of their names in literatures, in order to facilitate clinical application of these prescriptions (i. e. unification of their Arabic and Chinese names, and textual research of identical drugs with different Arabic names). This lays a foundation for the clinical application of Hui prescriptions and the analysis on compatibility regulatory, and provides scientific basis for studies on new Hui medicines.

Multiparametric bifurcations of an epidemiological model with strong Allee effect.

In this paper we completely study bifurcations of an epidemic model with five parameters introduced by Hilker et al. (Am Nat 173:72-88, 2009), which describes the joint interplay of a strong Allee effect and infectious diseases in a single population. Existence of multiple positive equilibria and all kinds of bifurcation are examined as well as related dynamical behavior. It is shown that the model undergoes a series of bifurcations such as saddle-node bifurcation, pitchfork bifurcation, Bogdanov-Takens bifurcation, degenerate Hopf bifurcation of codimension two and degenerate elliptic type Bogdanov-Takens bifurcation of codimension three. Respective bifurcation surfaces in five-dimensional parameter spaces and related dynamical behavior are obtained. These theoretical conclusions confirm their numerical simulations and conjectures by Hilker et al., and reveal some new bifurcation phenomena which are not observed in Hilker et al. (Am Nat 173:72-88, 2009). The rich and complicated dynamics exhibit that the model is very sensitive to parameter perturbations, which has important implications for disease control of endangered species.

[Experimental research of the inhibition of vascular endothelial growth factor C expression in gastric cancer by targeting RNA interference].

OBJECTIVE: To construct the plasmid expression vector pSIH1-H1-copGFP for RNA interference against vascular endothelial growth factor C (VEGF-C) and to evaluate its effect on the expression of VEGF-C mRNA in gastric cancer cells after transfection. METHODS: Three siRNAs of genome sequence of VEGF-C gene were retrieved from GenBank and one negative chain was used as control. Four siRNAs were cloned into plasmid pSIH1-H1-copGFP,which were then transfected into gastric cancer cells (SGC7901). The expression of VEGF-C mRNA was analyzed by RT-PCR. RESULTS: The recombinant plasmid of pSIH1-H1-copGFP specific for VEGF-C was confirmed by gene sequencing analysis. The target sequence obtained was completely consistent with the design. Transfection efficiency of the three siRNAs ranged from 60% to 70%. After transfection, the expression of VEGF-C mRNA in SGC7901 cells was significantly inhibited. Inhibition rates of VEGF-C mRNA expression were 35.4%, 33.8% and 81.5% in the three siRNA plasmid vectors, respectively. CONCLUSION: The siRNA expression plasmid vector against VEGF-C mRNA is successfully constructed, and RNAi may be a useful technique to inhibit the lymphangiogenesis of gastric cancer.